Angular variation of hard back-to-back hadron suppression in heavy-ion collisions 
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The basic idea of jet tomography is to infer information about the density evolution of the medium 
created in heavy-ion (A-A) collisions by studying the suppression of hard probes in an A-A environ- 
ment as compared to the baseline process known from p-p collisions. The suppression of back-to- 
back correlations in heavy-ion collisions allows, due to a different geometrical bias, a view into the 
medium which is qualitatively different from the one offered by single hadron suppression. A control 
parameter for the suppression corresponding to a systematic variation of in-medium pathlengths and 
density can be obtained by studying collisions at finite impact parameter 6. A systematic variation 
of pathlength can then be introduced by studying the suppression pattern as a function of the an- 
gle cj> with the reaction plane. Using a 3-d hydrodynamical evolution model for the medium and a 
Monte-Carlo (MC) model which has been shown to successfully reproduce the measured suppression 
of back-to-back correlations in central collisions of Au-Au at 200 AGeV, we compute the suppression 
as a function of <j> for 6 of 2.4 fm. 4.5 fm, 6.3 fm and 7.5 fm. Given that this involves variations 
in both control parameters b and 4> a comparison with data should eventually allow to place strong 
constraints on the combination of energy loss model and medium evolution model. 

PACS numbers: 25.75.-q,25.75.Gz 



I. INTRODUCTION 

The original hope of studying medium-induced modifica- 
tions to hard processes in heavy-ion collisions has been 
to do jet tomography on the medium [l|, [3, H, 0, H, @- 
The underlying idea is as follows: Hard processes take 
place simultaneously with soft processes responsible for 
the creation of soft bulk matter. This places partons 
emerging from hard vertices inside the medium, and in- 
teractions with the soft medium subsequently lead to an 
energy transfer from hard partons to the soft medium, 
resulting in a suppression of observed hard hadron yield. 
The strength of the interaction with the medium reflects 
the density of the medium, thus one should be able to 
gain information about the medium density from the ob- 
served strength of the suppression. The chief observable 
considered so far has been the nuclear suppression factor 
of single hadrons Raa (cf. e.g. @, Q) which is the mea- 
sured yield in A-A collisions divided by the yield in p-p 
collisions multiplied with the number of binary scatter- 
ings (i.e. the default expectation if there would not be 
soft processes forming a medium in an A-A collision) . 

Unfortunately, there are a number of problems when one 
tries to carry out this program in practice. First, Raa is 
not a sensitive measure of the averaged energy loss proba- 
bility distribution (P(AE, E)) Taa Assuming that 
eikonal propagation is a good approximation for hard 
partons, this quantity is the complete momentum space 
information about energy loss — given that the vertex of 
origin of any observed hard hadron cannot be known, pro- 
duction point, direction and hence the length of propaga- 
tion through the medium must be averaged probabilisti- 
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cally, leaving a momentum space probability distribution 
that a parton with original energy E lost the amount AE 
to the medium. The fact that (P(AE, E))t aa cannot be 
reliably determined from Raa is rather unfortunate, as 
knowledge of the energy loss probability distribution in 
momentum space is necessary before any mapping into 
position space (and hence any determination of density 
profiles) can be attempted. 

The second obstacle is that there is not even qualita- 
tive agreement what the precise nature of in-medium 
energy loss is. In addition to medium-induced radia- 
tion P, 0, H, El IE [fit, a number of publications have 
advocated a sizeable component of elastic energy loss 
EE 03 0, Q Ell which has a qualitatively dif- 
ferent pathlength dependence (linear in a homogeneous 
medium, as opposed to quadratic for radiative energy 
loss). As energy loss models typically contain one ad- 
justable parameter linking the medium density with the 
interaction strength, Raa could be reproduced by a num- 
ber of very different scenarios. 

Some steps can be taken to improve the situation. First, 
instead of trying to determine an a priori unknown den- 
sity evolution using hard probes, one may settle for the 
more modest goal of using a medium evolution model 
which is already constrained by soft hadronic observables 
and study the variation of the hard probes with changes 
of a known control parameter such as centrality. In this 
way, one employs hard probes as a test of a specific com- 
bination of evolution and energy loss model. Using the 
same 3-d hydrodynamical evolution which has been suc- 
cessful in describing bulk matter evolution [l7|, it has 
been shown that different energy loss models predict a 
different suppression pattern as a function of reaction 
plane when going to off-central collisions even when all 
are adjusted to describe the observed suppression for cen- 
tral collisions 0, US, H(J . 
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One may also turn to observables in which the medium 
geometry is probed in a qualitatively different way, such 
as hard back-to-back correlations. Such correlations have 
been measured by the STAR collaboration [21] . Employ- 
ing again the strategy to use a set of a priori constrained 
medium evolution models, it could be shown in [22I I23I ] 
that the pathlcngth dependence of radiative energy loss is 
compatible with the strength of back-to-back suppression 
whereas a linear pathlength dependence characteristic of 
clastic energy loss is not [13, thus any contribution with 
linear pathlcngth dependence must be a small correation. 
In the present paper, we continue the work of [H, H3 in 
computing the suppression of back-to-back correlations 
for non-central collisions as a function of the angle with 
the reaction plane. Since this is a very differential ob- 
servable, it constitutes a rather sensitive test of the ra- 
diative energy loss model 0, 0] in combination with the 
3d-hydrodynamical evolution model . 



II. THE MODEL 

The main ingredients of the model we wish to use have 
been extensively described in [13, [H, [H, [H, [25]. We 
will therefore confine ourselves here to a summary and 
describe only the improvements made to the framework 
in more detail. 

We calculate the correlation strength of hadrons back to 
back with a hard trigger in a MC simulation. There are 
three important building blocks to this computation: 1) 
the primary hard parton production, 2) the propagation 
of the partons through the medium and 3) the hadroniza- 
tion of the primary partons. Only step 2) probes prop- 
erties of the medium, and hence it is here that we must 
specify details of the evolution of the medium and of the 
parton-medium interaction. 

The strength of the parton-medium interaction contains 
one adjustable parameter K (see below). This parame- 
ter is fixed by the requirement that the model should de- 
scribe the suppression of single inclusive pions for central 
collisions. The results for non-central collisions, back-to- 
back correlations, different hadron species and different 
orientation with respect to the reaction plane are (given 
the hydrodynamical calculation as input) obtained with- 
out additional free parameters. 



A. Primary parton production 

The production of two hard partons fc, I in leading order 
(LO) perturbative Quantum Choromdynamics (pQCD) 
is described by 

d(J AB^kl + X ^^ kl 

dp T dy x dy 2 ~" 1 dt 

(i) 



where A and B stand for the colliding objects (protons 
or nuclei) and 2/1(2) is the rapidity of parton k(l). The 
distribution function of a parton type i in A at a mo- 
mentum fraction X\ and a factorization scale Q ~ px is 
fi/A{x\,Q 2 ). The distribution functions are different for 
the free protons [2(| [2tJ and nucleons in nuclei [H Hi| . 
The fractional momenta of the colliding partons i, j are 
given by x 1<2 = ^ (exp[±yi] + exp[±y 2 ])- 

Expressions for the pQCD subprocesses ^-ji — (s, t, u) as 
a function of the parton Mandelstam variables s, t and 
u can be found e.g. in [3(|. By selecting pairs of k, I 
while summing over all allowed combinations of i,j, i.e. 
99i 91-, 9Qi 99i <7<Z> 99 where q stands for any of the quark 
flavours it, d, s we find the relative strength of different 
combinations of outgoing partons as a function of px- 
For the present investigation, we require 2/1=2/2=0, 
i.e. we consider only back-to-back correlations detected 
at midrapidity. In a first step, we sample Eq. {1} summed 
over all k,l to generate px for the event, in the second 
step we perform a MC sampling of the decomposition of 
Eq. (TTJ) according all possible combinations of outgoing 
partons k, I at the px obtained in the first step. We thus 
end with a back-to-back parton pair with known parton 
types and flavours at transverse momentum px ■ 
To account for various effects, including higher order 
pQCD radiation, transverse motion of partons in the nu- 
cleon (nuclear) wave function and effectively also the fact 
that hadronization is not a collinear process, we fold into 
the distribution an intrinsic transverse momentum kx 
with a Gaussian distribution, thus creating a momentum 
imbalance between the two partons as pxi + PT2 = ^t- 



B. Parton propagation through the medium 

The probability density P(xo,yo) for finding a hard ver- 
tex at the transverse position ro = (xo,yo) and impact 
parameter b is in leading order given by the product of 
the nuclear profile functions as 

p( s T A (r + h/2)T A (v -h/2) 

p( ^ o) = iu4b) ' (2) 

where the thickness function is given in terms of 
Woods-Saxon the nuclear density p^(r, 2) as Ta(i-) = 
/ dzpA{r, z). Note that Eq. © may receive (presumably) 
small corrections when going beyond a leading order cal- 
culation. 

If we call the angle between outgoing parton and the 
reaction plane <f>, the path of a given parton through the 
medium £(t) is specified by (ro, (j)) and we can compute 
the energy loss probability P(AE) pat h for this path. We 
do this by evaluating the line integrals 

/>oo />oo 

u c {r ,4>)= dtiq(0 and (qL)(r , <j>) = / d&(0 
Jo Jo 

(3) 
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along the path where we assume the relation 

= K ■ 2 • e 3/4 (£) (cosh p - sinh p cos a) (4) 

between the local transport coefficient q(£) (specifying 
the quenching power of the medium) , the energy density 
e and the local flow rapidity p with angle a between flow 
and parton trajectory [3Ul32|. Energy density e and local 
flow rapidity p are the input from the 3-d hydrodynamical 
simulation of the medium evolution [r?| . 
u> c is the characteristic gluon frequency, setting the scale 
of the energy loss probability distribution, and (qL) is a 
measure of the path-length weighted by the local quench- 
ing power. We view the parameter K as a tool to account 
for the uncertainty in the selection of a s and possible 
non-perturbative effects increasing the quenching power 
of the medium (see discussion in [23|]) and adjust it such 
that pionic Raa for central Au-Au collisions is described. 
This leads to a value of K = 3.6 

Using the numerical results of [331 ] . we obtain 
P(AE; uj c , R) pa th for ^ c and R = 2uj 2 c /{qL) for given jet 
production vertex and angle <f>. In the MC simulation, 
we first sample Eq. ^ to determine the vertex of origin. 
For a given choice of 4>, we then propagate both partons 
through the medium evaluating Eqs. ([3]) and use the out- 
put to determine P(A£';w c , R) pa th which we sample to 
determine the actual energy loss of both partons in the 
event. 



C. Hadronization 

Finally, we convert the simulated partons into hadrons. 
More precisely, in order to determine if there is a trig- 
ger hadron above a given threshold, given a parton k 
with momentum pt, we need to sample A\~^ h {zi,pT), 
i.e. the probability distribution to find a hadron h from 
the parton k where h is the most energetic hadron of the 
shower and carries the momentum Pt = z\ -px- In the 
following, we make the assumption that the hadroniza- 
tion process itself, at least for the leading hadrons of a 
shower, happens well outside the medium and as a con- 
sequence neglect any interaction of formed hadrons with 
the medium. The timescale for hadronization of a hadron 
h in its restframe can be estimated by the inverse hadron 
mass, r/j ~ 1/nih, boosting this expression to the lab 
frame one finds Th ~ E^/mi. Inserting a hard scale of 6 
GeV or more for the hadron energy and the pion mass in 
the denominator (as pions constitute the bulk of hadron 
production), this assumption seems well justified. 
In previous works (22l . [23| we have approximated this 
by the normalized fragmentation function Dk-*h(z, Pt), 
sampled with a lower cutoff z m i n which is adjusted to the 
reference d-Au data. This procedure can be justified by 
noting that only one hadron with z > 0.5 can be pro- 
duced in a shower, thus above z = 0.5 the Dk^h{z, Pr) 
and Ai~* h (zi,pr) are (up to the scale evolution) iden- 
tical, and only in the region of low z where the frag- 



mentation function describes the production of multiple 
hadrons do they differ significantly. 

We improve on these results by extracting Ai(zi,pt) 
and the conditional probability to find the second most 
energetic hadron at momentum fraction zi given that 
the most energetic hadron was found with fraction z\ 
A-i{z\, Z2,Pt) from a shower evolution code. The proce- 
dure is described in detail in [25j where we used PYTHIA 
[HI to simulate the shower, in the present paper we em- 
ploy HERWIG [111 instead which provides a slightly bet- 
ter description of the d-Au baseline data. A comparison 
of both approaches with the d-Au data is shown in Fig. [I] 
Our way of modelling hadronization corresponds to an 
expansion of the shower development in terms of a tower 
of conditional probability denities Ajy (zi , . . . , z n , p) with 
the probability to produce n hadrons with momentum 
fractions Z\ , . . . z n from a parton with momentum px be- 
ing n™ =1 Ai(zi, . . . Zi,pr). Taking the first two terms of 
this expansion is justified as long as we are interested in 
sufficiently hard correlations. 

Sampling A\{z\,pt) for any parton which emerged with 
sufficient energy from the medium provides the energy 
of the two most energetic hadrons on both sides of the 
event. The harder of these two is by definition the trigger 
hadron and defines the near side. The hadron opposite 
to it is then the leading contribution to the away side 
correlation. 

The leading contribution to the correlation strength on 
the trigger side arises only in NL fragmentation, i.e. when 
sampling Ai(z\,zi,Pt) for the (fixed) z\ of the trigger 
to find the momentum of the second most energetic near 
side hadron. This term provides a correation to the away 
side correlation strength. This correction is Pt depen- 
dent and decreases from 25% in the lowest Pt bin in the 
present investigation (4-6 GeV) to 2.7% in the highest 
(10+ GeV) bin. Since we make a model comparison with 
the baseline both on near and away side, we compute 
both leading term and next-to-leading corection on the 
away side for consistency. 



D. Hydrodynamic Description of the Medium 

For the hydrodynamical description of the soft medium, 
we use the 3-d rclativistic fluid dynamics model by Bass 
and Nonaka [l7j]. The starting point is the relativistic 
hydrodynamic equation 

d^ v = 0, (5) 

where T MI/ is the energy momentum tensor which is given 
by 

T" u = (e + p)U»U v -pg» v . (6) 

Here e, p, U and g^ v are energy density, pressure, four 
velocity and metric tensor, respectively. The relativis- 
tic hydrodynamic equation Eq. ([5]) is solved numerically 
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FIG. 1: (Color online) Comparison of the 200 AGeV d-Au baseline measurement [21J for the near (left panel) and away 
side (right panel) correlation strength with calculations utilizing PYTHIA [34j or HERWIG [35| for shower development and 
hadronziation. 



using baryon number n B conservation 

d tl {n B (T^)U' M ) = Q. (7) 

as a constraint and closing the resulting set of partial 
differential equations by specifying an equation of state 
(EoS): e = e(p). In [l7| a bag model EoS is used for the 
QGP. 

The code utilizes a Lagrangian mesh and light-cone co- 
ordinates (r, x, y, rj) (r = \lt 2 — z 2 ), in order to optimize 
the model for ultra-relativistic regime of heavy collisions 
at RHIC. It is assumed that hydrodynamic expansion 
starts at tq = 0.6 fm. Initial energy density and baryon 
number density are parametrized by 

nB(x,y,v) = n Bmax W(x,y;b)H(r]), (8) 

where b and e max (^Bmax) are the impact parame- 
ter and the maximum value of energy density (baryon 
number density), respectively. W(x,y;b) is given 
by a combination of wounded nuclear model and bi- 
nary collision model and H(r]) is given by H(-q) = 
exp [— (|??| — ?7n) 2 /2g^ ■ 9(\ri\ — Vn)] ■ All parameters of 
the hydrodynamic evolution [ljj have been fixed by a 
fit to the soft sector (elliptic flow, pseudo-rapidity dis- 
tributions and low-pr single particle spectra), therefore 
providing a fully determined medium evolution for hard 
probes to propagate through. 

Note that the fluid description in principle contains also 
hard hadrons from the high Pt tails of the thermal dis- 
tributions. However, these contributions are small, and 
while the soft medium may influence hadron production 
between 4 and 6 GeV due to recombiantion processes, 
hadron production above 6 GeV is clearly dominated 
by hard processes [H, [3?], [H[ . Consequently we neglect 
mechanisms of hadron production other than hard pro- 
cesses in the following. 



III. RESULTS 

We compute the suppression of the away side correlation 
strength for a momentum range of 12-20 GeV for the 
trig ger hadron. This choice is motivated by the results of 
|23j where we found that information about the medium 
density becomes apparent only if the spread between trig- 
ger momentum and lowest associate momentum bin is 
large enough to allow not only particles with no energy 
loss to be registered but also include contributions from 
partons after finite energy loss. On the other hand, we 
found that the assumption that high Pt hadron produc- 
tion can be modelled by jet fragmentation only holds for 6 
GeV and above. Below, coalescence processes [361 l37l l38j 
play a significant role in hadron production in heavy-ion 
collisions. 

This requirement essentially determines the choice of the 
associate momentum windows: We include the window 
from 4-6 GeV for comparison, stressing that we expect 
that the model underpredicts the yield in this region. 
The remaining bins are chosen as 6-8 GeV, 8-10 GeV 
and 10+ GeV. In order to cancel systematic errors in the 
calculation of the baseline reaction without a medium, we 
present our results in the form of Iaa, i-e. the per-trigger 
yield in Au-Au collisions (for given b and <j)) divided by 
the per-trigger yield in p-p collisions. 
Since we do not find any significant modifications on the 
near side (cf. the discussion in [231]). we confine ourselves 
in the following to studying the modifications seen on the 
away side. 



A. Iaa as a function of Pt, b and (f> 

We show the resulting Iaa as a function of associate 
hadron momentum bin in Fig. [2] for four different val- 
ues of impact parameter b both for in-plane (0 = 0) and 
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FIG. 2: (Color online) Pt dependence of charged hadron Iaa 
in 200 AGeV Au-Au collisions for a 12-20 GeV trigger as a 
function of impact parameter b, shown for in-plane and out- 
of-plane emission. 

The first feature observed is that there is no strong Pt 
dependence observed in Iaa for any impact parameter, 
albeit a small rise may be present. The strongest trend 
seen is that Iaa shows a rise with b, corresponding to the 
fact that both entropy production (and hence medium 
density) and average pathlength are reduced for periph- 
eral collisions as compared to central collisions. A grow- 
ing split between in-plane and out-of-plane emission as 
a function of b shows that the spatial asymmetry in the 
initial state of the hydrodynamical calculation is indeed 
translated into a <fi dependent suppression of the correla- 
tion strength. 
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FIG. 3: (Color online) Pt dependence of the spread between 
in-plane and out-of-plane emission in 200 AGeV Au-Au colli- 
sions for a 12-20 GeV trigger as a function of impact param- 
eter b. 

We show the Pt dependence of the difference between 



in-plane and out-of-plane emission in Fig. [3] There is 
some indication that the separation between the differ- 
ent angular orientations is more pronounced at lower Pt, 
although the significant reduction of MC statistics in the 
high Pt bins as compared to the low Pt bins poses dif- 
ficulties to a solid statement. 

Surprisingly, the magnitude of both Iaa and of the split 
between in-plane and out-of-plane emission as a function 
of b is rather similar to the values found for single hadron 
suppression in terms of Raa found within the same for- 
malism [l8j. However, Raa and Iaa are in fact not the 
same, the similarity holds only on the level of about 30%. 
For this particular set of trigger and assoicate momenta, 
the magnitude of Iaa is larger than Raa whereas the 
spread between in-plane and out-of-plane emission is con- 
sistently smaller than the spread found in Raa- 



B. The geometry of single hadron and dihadron 
suppression 

We show the probability density of vertices in the (x, y) 
plane leading to a near side trigger hadron between 12 
and 20 GeV in Fig. 2] for impact parameters b = 2.4 
fm and b = 7.5 fm both for in-plane and out-of-plane 
emission. 

For the most central impact parameter b — 2.4 fm 
the results of the simulation in the 3-d hydrodynamical 
medium evolution are very similar to the results obtained 
in [23| for central Au-Au collisions in a 2-d hydrodynam- 
ics simulation (39[. There is a clear surface effect visible, 
the most likely point of origin of a triggered hadron is not 
in the medium center but about 4 fm displaced and hence 
close to the surface. The detailed position of this max- 
imum arises from a balance between the distribution of 
primary vertices coming from the nuclear overlap Eq. ([2|) 
which peaks in the medium center and the fact that par- 
tons originating in a high-density region are more likely 
to be quenched. On the other hand, close to the medium 
edge where the quenching is small, the probability to find 
a vertex is much reduced. 

For b — 2.4 fm, there is no large difference between in- 
plane and out-of-plane visible, corresponding to the fact 
that the spatial asymmetry of the medium is small. This 
is different for the simulation with b — 7.5 fm. Here, clear 
changes in geometry when going from in-plane to out-of 
plane emission are visible. In particular, the surface effect 
is stronger for out-of-plane emission due to the greater 
amount of medium in the path of outgoing partons for 
this configuration. 

Note that the distribution shown here is qualitatively (i.e. 
up to small changes with Pt) the same for any hard sin- 
gle hadron observable, in particular Raa (averaged over 
the trigger range) and 7-hadron correlations (where the 
photon escapes into the +x direction). 
We show the probability density of vertices in the (x, y) 
plane leading to a near side trigger hadron between 12 
and 20 GeV and an associate away side hadron with 4 
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b = 2.4 fm, in plane b = 2.4 fm, out of plane 
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FIG. 4: Probability density for an inclusive hard hadron measurement, i.e. for finding a parton production vertex at (x, y) given 
a triggered event with 12 GeV < pr < 20 GeV for different b and <f>. In all cases the near side (triggered ) hadron propagates 
to the —x direction. Countours are at linear intervals. 



GeV < px < 6 GeV in Fig. [5] It is immediately obvious 
that the distribution is very different from the distribu- 
tion of vertices for single hadron observation shown in 
Fig. 2J First, the dihadron distributions are much wider 
in ±y direction, indicating the importance of the periph- 
ery where both near and away side parton have a short 
in-medium path (or the halo where the production ver- 
tex lies outside the medium). Second, the distribution 
is almost symmetric with respect to the y-axis, indicat- 
ing the events are favoured in which the near and away 
side pathlength is about equal. Both features are consis- 
tent with the distributions observed in [23[ for different 



density evolution models. 

In particular, note the absence of tangential emission, 
which in the framework presented here does occur, albeit 
at much stronger quenching power of the medium (23j . 

It is evident from the underlying geometry that the 
physics (in terms of average pathlength and relevant 
medium density) of single hadron and dihadron suppres- 
sion is rather different in the calculation. Geometry alone 
clearly cannot account for the observed similarity be- 
tween Raa and Iaa- 
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b = 2.4 fm, in plane b = 2.4 fm, out of plane 
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FIG. 5: Probability density for a coincidence hard hadron measurement, i.e. for finding a vertex at (x,y) leading to a triggered 
event with 12 GeV < pr < 20 GeV and in addition an away side hadron with 4 GeV < pr < 6 GeV for different b and cj>. In 
all cases the near side hadron propagates to the —x direction. Countours are at linear intervals. 



IV. DISCUSSION 

One may be tempted to conclude from the similarity of 
Raa and Iaa that there is no additional information con- 
tained in Iaa beyond what can be obtained from single 
hadron suppression. However, this is clearly not the case. 
Let us illustrate in a qualitative way the meaning of this 
observation. 

In a purely geometrical suppression picture, a parton is 
absorbed when the hard vertex of its origin is in a dense 
region but unmodified when the vertex falls into the di- 



lute halo region. In this case, the whole back-to-back 
event is either absorbed or unmodified. This in turn im- 
plies that the number of triggered events is reduced, but 
once an event is triggered, there is no modification of 
the near or away side parton. From this one may imme- 
diately deduce that in this scenario Raa may be below 
1, but Iaa will always be (modulo small effects by the 
nuclear parton distribution functions) equal to 1. 
Thus, in any situation in which quenching is chiefly 
driven by the density at the vertex of origin and the 
medium can be modelled as an opaque and a dilute re- 
gion, Iaa is expected to be substantially larger than 
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Raa- Within this model, this trend has been demon- 
strated as a rise in the per-trigger yield in the black core 
scenario described in [231 ]. 

Let us next consider the momentum space. By definition, 
the trigger hadron is the hardest hadron of the event. 
Consequently, there is the strongest bias on the trigger 
parton to experience no or only a small amount of en- 
ergy loss. In contrast, the bias given that there is a hard 
hadron on the near side on the away side hadron is con- 
siderable reduced, allowin g av erage energy losses of sev- 
eral GeV on the away side [4(| . If the geometry-averaged 
energy loss probability on near and away side would be 
equal (i.e. if there would not be a bias towards surface 
emission), this would imply again that Iaa should be 
substantially larger than Raa (for the measured Raa of 
about 0.2, this would amount to almost a factor 2 differ- 
ence). 

However, there is an opposing bias in position space: The 
fact that the hardest hadron tends to be produced close 
to the surface means that the away side parton on average 
has a longer pathlength. Here, the weighting of the en- 
ergy loss with pathlength is crucial, as the medium den- 
sity drops due to the longitudinal and radial expansion of 
the medium while the parton traverses the medium. An 
energy loss which grows linear in pathlength L in a con- 
stant medium (as characteristic of elastic energy loss) for 
example translates into a logarithmic pathlength depen- 
dence ~ \xl(L/Lq) in a medium undergoing longitudinal 
Bjorken expansion where the density drops like 1/r. In 
this situation, pathlength differences on near and away 
side do not matter much and the bias in momentum space 
is expected to win out. This is observed within the model 
described here in [24| . 

On the other hand, for an energy loss with quadratic 
pathlength dependence ~ L 2 in a constant medium, the 
energy loss gets proportional to ~ L in a medium under- 
going longitudinal Bjorken expansion. If the away side 
path is on average twice as long as the near side path this 
effect would be expected to cancel the reduced bias in mo- 
mentum space and make Raa and Iaa about equal. This 
argument strongly relies on modelling the expansion of 
the medium correctly — note that in the static medium 
studied in [41] in a conceptually similar framework us- 
ing the same radiative energy loss pciture, back-to-back 
events are almost completely suppressed, corresponding 
to an energy loss dependence with pathlength which is 
too strong. 

These arguments may help to understand why the fact 
that Iaa and Raa are of similar magnitude is in fact 
highly non-trivial and in itself probes properties of the 
medium and its expansion along with properties of the 
parton-medium interaction. 



V. CONCLUSIONS 

We have computed the variation of the suppression of 
back-to-back correlations for non-central 200 AGeV Au- 
Au collisions as a function of impact parameter b for both 
in-plane and out-of-plane emission. Since the parame- 
ter determining the quenching power of the medium has 
been adjusted such that single hadron Raa for central 
collisions is described, the extension to back-to-back cor- 
relations, non-central collisions and the dependence on 
the angle with the reaction plane is computed without 
additional free parameters, given the underlying hydro- 
dynamical evolution (which is in turn constrained by the 
need to describe bulk matter observables such as mo- 
mentum spectra and elliptic flow v 2 ). Thus, the results 
represent a rather constrained prediction which can be 
used to test the validity of the combination of hydro- 
dynamical evolution [1 71 ] with the radiative energy loss 
model @, HI. 

The suppression pattern we find is not unexpected given 
previous results [HI, [23[: As a function of associate 
hadron momentum, Iaa is approximately flat but may 
exhibit a small rise. Iaa increases with impact param- 
eter, reflecting the reduced soft matter entropy produc- 
tion and a shortening of average in-medium pathlength. 
At the same time, a split between in-plane and out-of- 
plane emission grows, reflecting the spatial asymmetry 
in the initial state. There is an interesting similarity be- 
tween the numerical values of Iaa (and the magnitude 
of the split) and the values of Raa and the split be- 
tween in-plane and out-of-plane emission observed there. 
There is no single physics reason which would dictate 
this outcome, rather one is looking at a cancellation of 
various different effects. An investigation in the geome- 
try reveals that the physics underlying single hadron and 
dihadron observables in the model is rather different in 
terms of dominant point of origin, typical density being 
probed and average pathlength. A detailed comparison 
with data and computations for different system size, e.g. 
Cu-Cu may be helpful to investigate this further. 
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